Spin chains of Haldane-Shastry type and a generalized central limit theorem 
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We show that the density of energy levels of a wide class of finite-dimensional quantum systems 
tends to a Gaussian distribution as the number of degrees of freedom increases. Our result is 
based on a nontrivial modification of the classical central limit theorem, and is especially suited to 
models whose partition function is explicitly known. In particular, we provide the first theoretical 
explanation of the fact that the level density of several spin chains of Haldane-Shastry type is 
asymptotically Gaussian when the number of sites tends to infinity. 
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Spin chains of Haldane-Shastry (HS) type are the 
prime example of integrable spin chains with long-range 
interactions, having close connections with several topics 
of current interest such as strongly correlated systems , 
generalized exclusion statistics 4, 5j, and the AdS-CFT 
correspondence [6] . Recent numerical computations have 
shown that a common feature of these chains HJl, S 
and their supersymmetric extensions 11 , [l^. Il3l| is the 
fact that the level density (normalized to 1) becomes 
Gaussian when the number of sites N tends to infinity. 
The knowledge of the continuous part of the level density 
is a key ingredient in the theory of quantum chaos, since 
it is used to rescale ( "unfold" ) the spectrum as a prelimi- 
nary step in the study of important statistical properties 
such as the distribution of spacings between consecutive 
levels 1J|. In fact, a long-standing conjecture of Berry 
and Tabor [l^ posits that the spacings distribution of a 
"generic" integrable system should be Poissonian, while 
in quantum chaotic systems like polygonal billiards this 
distribution is given by Wigner's surmise flG^, character- 
istic of the Gaussian ensembles in random matrix the- 
ory. For many spin chains of HS type, it can be shown 
that the Gaussian character of the level density implies 
that the spacings distribution obeys neither Poisson's nor 
Wigner's law, but is rather given by a sirnple "square- 
root of a logarithm" formula [lOi ^^^^ 
letter we develop a generalization of the standard cen- 
tral limit theorem to show that the level density of a 
wide class of finite-dimensional quantum systems must 
be asymptotically Gaussian. This class includes in par- 
ticular a supersymmetric version of the original (trigono- 
metric) Haldane-Shastry spin chain, as well the BCn 
version of the spin 1/2 Polychronakos-Frahm (ratio- 
nal) chain [13, [ll- 



To explain the gist of our approach, let us begin by 
considering an extremely simple toy model for which 
the asymptotically Gaussian character of the spectrum 
can be easily established. Indeed, let H he a finite- 
dimensional Hilbert space and let us denote by 7i®^ its 



N-th tensor product, which can be identified with the 
TV-particle space. Our toy Hamiltonian will be given by 

H = h + --- + lN, (1) 

where each operator Ik = 1^^ • • • (g) / (g) • • • (g) 1-h is a con- 
stant of motion acting on the Hilbert space of the A:-th 
particle as some fixed operator /. Physically, this sys- 
tem can be thought of as consisting of N identical, non- 
interacting particles. Had we allowed the one-particle 
Hilbert space to be infinite dimensional, this would pre- 
cisely be the structure underlying the Ai"-dimensional har- 
monic oscillator. 

It is not difficult to see that in the limit N oo the 
level density of the Hamiltonian ([1]) becomes Gaussian, 
with mean /i = N^i and variance = Na'j {hj and 
cr| respectively denoting the mean and variance of the 



spectrum of /). 
by 



Indeed, the eigenvalues of H are given 



E 



N ■ 



(2) 



where each £k is an eigenvalue of the operator /. Since 
each of these eigenvalues can be freely chosen, the pa- 
rameters £k in the above formula can be regarded as in- 
dependent random variables with the same distribution. 
Thus, the fact that the random variable £ asymptotically 
follows the Gaussian law as — > oo is an immediate con- 
sequence of the central limit theorem. 

It is natural to wonder if the simple argument above 
can be extended to a wider class of models. A cursory in- 
spection reveals, however, that several assumptions must 
be substantially relaxed in order to cover any physically 
interesting situation. In particular, in Eq. ([2|) it is cru- 
cial to allow for sums of independent random variables 
which are not identically distributed. It is also clear that 
it would be preferable to express the conditions on the 
spectrum of our model in terms of its partition function 
rather than its eigenvalues. Indeed, for chains of HS type 
the partition function can be computed in closed form 
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through Polychronakos's "freezing trick" 17|, ll9|, while 
the spectrum is considerably more difficult to handle. 
Thus, in order to further develop the former approach, we 
will need to prove a modification of the classical central 
limit theorem incorporating the above requirements. 

Let us introduce some notation before going on. We 
shall denote by 



1=1 

the partition function of a finite-dimensional Hamiltonian 
H with energy levels £i and degeneracies dj, 1 ^ i ^ n, 
where q = e~^^^'^^^^ and fee is Boltzmann's constant. 
The value Z{1) = d-i is thus the dimension of the 
Hilbert space of the system. Its thermodynamic energy 
at temperature T is given by 

E{q)^q^\ogZ{q), 

in terms of which the mean and variance of the spectrum 
of H can be expressed as 

n 

^ ' 1=1 

We will henceforth restrict our attention to systems 
whose partition function can be written as a product 



N 



Z{q) = n Zk{q) , 



(3) 



fc=l 



and we shall assume that the factors Z^. possibly depend 
on the integer parameter N in view of forthcoming appli- 
cations. Roughly speaking, the system is equivalent to an 
effective model of N non-interacting, but not necessarily 
identical, subsystems. Two powerful tools in the study 
of the energy distribution of ([3]) are the characteristic 
function 



Z(e») 



(4) 



which is essentially its Fourier transform, and the nor- 
malized characteristic function 



(5) 



The key property of the latter function is that, under very 
mild technical assumptions [ISl , in the limit N ^ oo the 
energy distribution of the system asymptotically follows 
the Gaussian law with mean /i and standard deviation a 
if and only if 



N 



lim ip{t) 



(6) 



We shall next provide some simple conditions on the 
factors Zk ensuring that Eq. ([6]) holds. In order to state 
them in a concise way, we will use a subscript k (as in 
or (fk) to denote the statistical quantities computed with 
respect to the fc-th partition function Zk, and introduce 
the notation 



Mkir) 



sup 

|s|<r 



\0gipk{s) 



It should be noticed that the average and standard de- 
viation of the energy of the whole system are related to 
the analogous quantities of its subsystems by 



k 



(7) 



where the index runs from 1 to (as always hereafter, 
unless otherwise stated) . Our assumptions on the factors 
are then the following: 

(i) The standard deviation of the full partition func- 
tion and of its k-ih factor satisfy 



ICkCN, 



(8) 



where the constant Ci does not depend on N . 



(ii) There exist some positive constants ei, 62, C2, inde- 
pendent of N , such that 



Mk{ei) ^ CsiVs- 



(9) 



Let us briefly discuss the meaning of these assumptions. 
Roughly speaking, the first condition holds whenever the 
standard deviation of the whole system does not effec- 
tively depend only on a few subsystems. More precisely, 
let 

M{N) = max {cr^ : fc = 1, . . . , iV} , 

and suppose that there exists a constant a > (indepen- 
dent of N) and a function m{N) > such that 

al ^ miN) 

for at least aN values of k. Then, a sufficient condition 
ensuring the validity of (i) is that there exists a constant 
A such that 



M{N) ^ Am{N) 



(10) 



as N ^ 00. Indeed, by definition of m{N) and Eq. ([7]) 
we have that ^ aN ■m{N), which in turn implies that 

^ /ji£(iv) y . 

a \aNm{N) J ^ ^ ^ 
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with Ci = {A/aY/"^. On the other hand, the technical 
condition (ii) is a local bound on the third central mo- 
ment of the thermodynamic energy of the fc-th subsystem 
at imaginary temperature. Indeed, it can be shown that 



3 



1 d 



logZfc(e"fc) 



— {£^-Ekiq)y 



where the overbar denotes thermal average. 

It should also be noted that our second condition is 
reminiscent of Lyapunov's classical condition [23|, ac- 
cording to which ^ holds if there is a positive 5 such 
that 



lim 

N~*oo a 



(11) 



Note, however, that due to the absolute value the latter 
condition cannot be expressed in terms of the partition 
functions Zk (this is also true for the more general con- 
dition due to Lindeberg [13] )■ Thus pT|) is impractical 
when the partition function is known but there is not an 
explicit formula for the spectrum, as is the case with spin 
chains of HS type. 

In order to see that conditions (i) and (ii) above imply 
that the level density asymptotically follows the Gaussian 
law, let us compute the hmit of \og(p{t) as — > oo. To 
begin with, one should observe that the fact that (pk is 
normalized to zero mean and unit variance ensures that 
the second-order Taylor expansion of log (pk around is 



where the remainder is bounded by 



\RkiT)\ ^ -Mk{\r\). 



(12) 



Since (p{t) — Y[k Vk{t) by Eq. (jH), it immediately follows 
from © and ^ that 



iog^(0 = Ei°g^4v) = 4+E^4^)^ 



where by virtue of Eq. (fT^ the error can be estimated as 



log(^(i) + y 



(13) 



Let us assume that \t\ < to, where to is a fixed but other- 
wise arbitrary constant, and take N greater than (^^^) • 
In this case, by ^ we have 



<yk\t\ 



< CitoN-^ ^ ei. 



so that (IT3l) can be controlled as 



\og^{t) + - 



^ ^iCitofC2N-'' 
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on account of @. It then follows that \og(fi{t) converges 
pointwise to — for all real t, the convergence being 
uniform on compact sets. Hence Eq. ^ holds, and thus 
the level density becomes asymptotically Gaussian by the 
properties of characteristic functions, as we wanted to 
show. 

A particularly simple class of partition functions of the 
form ([3]) is obtained by requiring that each factor 
corresponds to a two-level system. In this case, setting 
the ground state energy to zero without loss of generality, 
each factor can be written as 



(14) 



and the mean and standard deviation of its energy are 
readily computed as fJ,k = cTk = \£{k,N). Its reduced 
characteristic function is simply Cpk{T) — cost, so that 
condition (ii) is automatically satisfied. It can be shown 
that condition (i) is satisfied as well, e.g., whenever 
£{k,N) depends polynomially on k and N. More pre- 
cisely, if r = deg£ then a'^ ^ N'^'^'^^, while is at most 

We shall next discuss how the previous developments 
can be directly applied to two specific spin chains of HS 
type, which are of considerable interest in themselves. 
The first one is the Polychronakos-Frahm chain of BCn 
type (2ll |. whose Hamiltonian is given by 



H 



E 



1 + eSjk 1 + eSjSkSjk 



(0-a-)' (0+6)^ 



e'Sk 



where e = e'~ = 1, /3 > 0, Sjk is the operator that 
permutes the j-th and k-th spins and Sk is the operator 
flipping the fc-th spin. Moreover, the chain site ^k is 
expressed in terms of the k-th zero yk of the generalized 
Laguerre polynomial as ^k — V^Uk- It has recently 

been shown [9| that for spin 1/2 the partition function of 
this model is given by 



Z{q) = q 



(15) 



where 6ie is the Kronecker delta. Other than the inessen- 
tial factor g2^(^^i)'5i= ^ which can be removed by shifting 
the ground state energy, this partition function is pre- 
cisely of the form (I14p with £{k, N) = k. It follows imme- 
diately from our previous discussion that when N oo 
the spectrum of H is normally distributed, with mean 
and variance given by 



E' 



N 



(N+l) 



E 



N 
12 



(7V+i)(7V+l) 
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This fact had been numerically verified in Ref. . 

The above result has an interesting interpretation in 
classical partition theory. Indeed, by Eq. (fTSj) the ener- 
gies of the ferromagnetic chain (e = —1) are the inte- 
gers in the range 0, 1, ... , N{N + l)/2, the degeneracy 
of an energy k being the number Qnik) of partitions of 
the integer k into distinct parts no larger than N (with 
(5jv(0) = 1). We have thus established the asymptotic 
formula 



QN{k) 



27rcr 



(with k = 0, 1, ... , N{N + l)/2, and ^, a given in the 
previous equation) which, to the best of our knowledge, 
was not previously known. 

The second spin chain we shall consider is the su- 
persymmetric version of the celebrated Haldane-Shastry 
chain, with Hamiltonian given by 11 1 



H = 



-y— 



ok 



j<k 



(16) 



where e = 1. The supersymmetric spin permutation op- 
erator Pjk acts on an element |si,...,SAr) of the spin 
basis as 



1 "^k ! 



. 



where ejfc(s) is —1 when either both Sj and Sk are 
fermionic spins, or Sj and Sk are spins of different type 
with an odd number of fermionic spins between them, 
and 1 otherwise. In the su(l|l) case (i.e., when there 
is only one bosonic and one fermionic internal degree of 
freedom), the partition function of the chain (fT6|) can be 
written as 12211 



z{q) ^ 2 n 



'1 + g'=(^-^-)) 



fe=l 



for both e = 1. Apart from the irrelevant factor of 2 
(which does not affect the level density, since it is nor- 
malized to 1), this partition function is a product of 
— 1 factors of the form (|14p with a polynomial func- 
tion £{k,N) = k{N ~ k). Hence the above argument 
rigorously establishes that in the limit N ^ oo the level 
density of the Hamiltonian (fTH]) becomes Gaussian, with 
parameters 



lj2k{N-k) 



N 
12 



1) 



a^ = \j:k^{N-kr^^^{N^-l) 



in whole agreement with the numerical computations in 
Ref. [iH. 



To conclude, let us summarize our results and offer 
some perspectives. We have provided conditions on the 
partition function of a finite-dimensional quantum sys- 
tem depending on a positive integer N ensuring that its 
level density is asymptotically Gaussian as TV ^ oo. Our 
conditions, which are related to Lyapunov's generaliza- 
tion of the classical central limit theorem, are directly 
formulated in terms of the partition function and do not 
require the explicit knowledge of the spectrum. We have 
applied our result to rigorously show that the level den- 
sity of two well-known spin chains of HS type becomes 
Gaussian as the number of sites tends to infinity. The 
first chain discussed is associated with the BCn root sys- 
tem and presents rational interactions, while the second 
one is a supersymmetric version of the original (Ajy-i- 
type) Haldane-Shastry spin chain. 

Our result does not apply to all spin chains of HS type, 
since the partition function of these models in general 
does not factorize as in Eq. What seems to be true 
in this case, however, is that there is an approximate 
factorization 



Z{q) = {l + e{N,q))l[Zk{q) 



(17) 



where the error term e(iV, q) and its g-derivatives tend 
to in a controlled way as the number of sites N tends 
to infinity. Although we shall not further elaborate on 
this point here, it is clear that the overall factor l-|-e(A^, q) 
can be taken into account in our approach by imposing 
appropriate technical conditions analogous to ([8]) and ([9]) . 
We shall provide a more detailed discussion of this issue 
in a forthcoming paper. 
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